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Abstract — A rigorous method for solving 2D scattering by an 
arbitrary perfectly conducting inclusion located in the region of 
a rectangular waveguide T-junction is presented. This method is 
developed from the approach based on the Green’s theorem. The 
latter is simultaneously the wave equation and the boundary 
conditions on the scatterer surface. The proposed method is 
illustrated by results obtained for a thick septum located inside 
the interaction region. Such approach can be applied for 
optimization of waveguide junctions.generalized using special 
weighting functions satisfying 

Index Terms — Waveguide Junctions, Green’s Theorem, 
Microwave, Rectangular Waveguide. 

I. INTRODUCTION 

The T-junction of two waveguides is an important element 
of waveguide tracts, which can be used as one of the key 
building blocks of various devices, e.g., power dividers, 
filters, multiplexers. The presence of extensive applications 
necessitates the development of methods for the study of such 
stmctures. 

The paper presents a rigorous solution for a 2D scattering 
problem of electromagnetic wave by a perfectly conducting 
inclusion of an arbitrary cross section shape within an 
interaction region of rectangular waveguide T-junction. In 
Fig. 1 the interaction region is bounded by line segments 
Lq , Ly , , and the scatterer surface L s . This model can 

be applied for optimizing the transmission properties of 
waveguide junctions by placing rods or septums. 

The approach proposed is based on the Green’s theorem 
method, which has been previously successfully applied to 
characterization of different T-junctions of rectangular 
waveguides. E.g. in [1] the problem is solved using the 
appropriate integral equation of the residue theory and the 
weighting functions, which were chosen so as to identically 
satisfy Helmholtz equation in the inhomogeneous region of 
the waveguide junction. The characteristic feature of our 
approach is that we choose the weighting functions in such 
manner that they automatically take into account the presence 
of the perfectly electric conducting (PEC) inclusion. More 
precisely, the system of wave functions is found from 
solutions of the three auxiliary problems for scattering by the 
inclusion in short-circuited waveguides. 
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Fig. 1. The structure under study: rectangular waveguide 
T-junction containing a PEC scatterer. 


A simple system of algebraic equations for determining the 
amplitudes of scattered fields in the waveguides can be 
obtained. To this aim only one propagating mode is taken into 
account in one of the waveguides and the scattering 
coefficients are found by using the three-short method [2, 3], 
which enables us to obtain the scattering coefficients with 
high accuracy. With this approach, one can construct the 
complete scattering matrix in a closed form for an arbitrary 
number of waveguide modes [4]. Thus, this paper presents the 
further development of the Green's theorem method for the 
case when the T-shaped interaction region of two variously 
sized rectangular waveguides is loaded with an E-plane 
perfectly conducting septum. 

II. FORMULATION OF THE PROBLEM 


For definiteness, we restrict our consideration to the 
analysis of LM waves in the E -plane of the waveguide 
junction [5]. The time dependence is assumed in the form 
exp {—ico • t) . For this type of waves the nonzero component 


E x (x 9 y 9 z) of the electric field can be presented in the 


following form: 

E x (x, y, z) = cos 


^ /mi ^ 
x 

V a j 


W(y,z). 


(l) 


The other projections of the electromagnetic field are written 
as follows: 
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where a is the common dimension of the waveguides. We 
assume that the X— mode order m = 0,1,2,... is fixed. 
Function W(y,z) satisfies the Helmholtz equation 

(A yz +K 2 )W(y,z)=0, (2) 

and zero boundary conditions on the PEC walls of the 
waveguide and contour L s of the scatterer S . Here, 
2 2 1/2 

K = [k - (nm / a) ] ' is the effective wavenumber, k 
is the wavenumber of the waveguide filling. In regular regions 
A(y>b ) , Z?(z<0) and C(z>b h ) , function 

W(y,z) is represented as 

+00 
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where and , C^ c are the amplitude 

coefficients of the incident and scattered waves; 
r/ n =[K 2 -(nn/b h ) 2 f 2 and 

/3 n = \K 2 — (nn / by ] 1/2 are the longitudinal 

wavenumbers of the respective waveguides. 

It is necessary to define the scattering matrix elements 
SU (I,J = A,B,C) , which are determined as 

A sc = SAA-A in + SAB • B in + SAC • C in , 

B sc = SBA A in + SBB B in + SBC C in , ( 3 ) 
C sc = SCAA in + SCB • B in + SCC • C in . 

Here we introduce the column vectors of the amplitudes of the 
incident and scattered waves A m = { } , B m = { B™ } , 

c in ={C;} ; A sc ={A; c } , B sc ={5-} , 

c sc ={c,n. 


Let’s introduce a solution to equation (3) in the 
inhomogeneous region W (x, y) that satisfies zero boundary 
conditions on the PEC surfaces: 


W 


= 0; W ^ =0, 


(4) 


Next we transform Helmholtz equation (3) inside the 
inhomogeneous region into the integral equation over the 

contour L = Lq+L [ +L 2 +L 3 +L s of the said region by 
applying the second Green’s formula [6]: 


~ dW dW 
dL\ W--W 


= 0, 


(5) 


dN dN 

V / 

where W(y,z) , as before, is an arbitrary solution to Eq. (3) 
in the inhomogeneous region and d / dN denotes the normal 


derivative at contour L . As W we consecutively substitute 
the three solutions of the auxiliary problems for the main and 
the side (see Fig. 1) waveguides. 

Taking into account the boundary conditions (5) on the 

perfectly conducting segments L 0 and L s of the integration 


contour the latter expression yields: 


dL 


L^+Lz+L^ 


w 


dW „dW 


dN 


W- 


dN 


= 0 . 


( 6 ) 


Next, we particularize our model and consider first the 
auxiliary problem, which represents the main waveguide 
containing the scatterer S in the form of a rectangular 
septum (see Fig. 2a). This model can formally be obtained 
from the structure presented in Fig. 1 if we consider the upper 

section of line Ly to be PEC. Let the scatterer in this auxiliary 
waveguide is excited by the mode of order q = 1,2,... 
propagating from z — — °o (port B). In this case, solution 
outside the scatterer W(y,z) =W^ (y,z) in region 


^3 < Z < Z 2 is 


W“(y,z) = 


WB n +W° c B ,z<z 2 , 

W£,z>z 3, 


(7) 


in equation (7) we have introduced the following notations for 
the waves: 


f nq ^ 

— y 

Vb y j 


W;=exp[^(z-z 2 )]sin 

+QQ 

^=LC ex PH/^-z 2 )]sin 

5=0 

+00 


, z<z 2 , 


f ns x 

AV 

f ns x 


s=0 


, Z<Z 2 , (8) 


Z > 


V u J 


As the second auxiliary problem, we consider the same 
main waveguide, but now assuming that the waveguide is 
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excited by the mode of order q = 1,2,... propagating from 
Z — +°° (port C). The solution outside the scatterer 
w(y,z) = w^(y, z ) looks like 


W 0 c (y,z) = 


I W?+W£,z>z 3 , 


W£,z<z 


(9) 


2 ’ 



(b) 


Fig. 2. The auxiliary structures: for the first and the second 
auxiliary problems (a); for the third auxiliary problem (b) 
where 
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Note that in (8), (10), and (12) S U (/, J = A, B, C ) are 
elements of the scattering matrices, which are defined during 
the auxiliary problems solution. 

Here we introduce the column vectors of the incident and 

scattered waves amplitudes: A in —{A^ 1 } , B in — {B ™} , 
C in ={C} ; A sc ={4* c } , B sc ={B s n c } , 


c sc ={ Cj 1 °}. The scattering coefficients and the solutions 

of the auxiliary problems can be found using various 
numerical techniques, such as the mode-matching method [7], 
boundary contour mode-matching method [8], the moment 
procedure [9, 10]. Here we use mode-matching method in the 
partial domains with PEC boundary conditions at the metal 

surfaces and continuity conditions for W(y,z) and its 

normal derivative either dW/dz or dW/dy on the relevant 
apertures. 

Next, we substitute the found weight functions W^(y, z) , 

W^(y, z) , Wq(y, z) to the second Green’s formula (6) in 

place of W and take into account the explicit form (4) for 
W(y, z) in each of the interaction regions and impose the 

boundary conditions for W and W on the PEC surfaces. As 
a result we arrive to a direct formula for calculating the 

scattering amplitudes A^ c , B^ c , C^ c : 


K +<■(«;; +«:)+<, -(C +oi, 

n - 1 

+00 

K = <+ + + ) + c i. <i3) 


n - 1 


+00 

<7 = 3 K ■ < ■ K + + ) + e”, ■ b: + oy • c i, 

n = 1 

where q = 1, 2, ... , and the matrix coefficients have the form: 


The third auxiliary problem is formulated for the structure 
formally obtained from Fig. 1 if the segments L 2 and L 3 are 
assumed to be PEC (see Fig. 2b). Let the wave of mode order 
q = 1,2,... propagates from the side y = +oo (port A). In 
this case the solution outside the scatterer 
W ( y, z) = Wq(y, z) for the region y > y 1 can be written 


R ;„ = exp[ " 7< )] n, exp [.77„ (b - y, )] • sy 


n. 


nq 


R-l = - — — GW^ B , 

qn 2 b h i qn 
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The integrals on the right hand sides of the above 
expressions can be calculated in a closed form provided that 

W q A ,W B ,W q are known. 

Next we substitute the second and third equations from (14) 

into the first equation thus eliminating {B^ c } and {C^ c }. 

As a result we obtain the system of linear algebraic equations 
(SLAE) for the amplitudes of the scattered field in the side 

waveguide A sc = { A^ c } , which in matrix form looks like 


(I - D) • A sc = U A • A in + IT • B in + IT • C in . (14) 

Here I - is identity matrix, and the other matrices appearing 
in (15) are defined as follows 


d = r b p a +r c q\ 

U a =R a +R b P a +R c Q a , 

U B =R B +R B P B +R C Q B , 

U c =R C +R B P c +R c Q c . 

The solution to SLAE (15) can be written in the form: 

A sc =(I-D) 1 -(U A -A in +U®B in +U c -C in ). (15) 


Comparing the solution (16) with the definition of the 
scattering matrix (13), we find expressions for the scattering 
matrix elements in the side waveguide A (y > b) 

S AA =(I-D)"‘-U\ 

S AB =(I-D) _1 -U b , (16) 

s AC = (i-d) _ 1 -u c . 


Using direct formulae (14-17) we obtain the elements of 
the scattering matrix in the main waveguide: 

• for scattering in port B of the main waveguide ( Z < 0 ) 

s BA =p A -(i-s AA ), 


gBB pH 


P A S 


AB 


(17) 


i^BC pC pA S AC 


• for scattering in port C of the main waveguide ( Z>b h ) 

S CA =Q A (I-S AA ), 

S CB =Q b -Q a S ab , (18) 

s cc =Q c -Q a S ac . 


III. NUMERICAL EXAMPLES 
Let us consider a PEC rectangular septum of thickness t 
and height h , placed in the interaction region of the empty 
waveguides and, to be specific, fixed to the downside wall of 
the main waveguide (see insert in Fig. 3). We are interested in 
the complex elements of the scattering matrix 

S U (/,/ = A,B,C) . The waveguide dimensions are 
a = 8.64cm , b = b h =a!2 . Figs. 3, 4 illustrate 
frequency dependence of the reflection coefficient for 
LM m impinged from the port A for various geometric 

parameters of the inclusion. In these figures the septum is 
placed symmetrically in the interaction region along the 
z-axis. 



Fig. 3. The reflection coefficient for LM 01 -mode at port A 
for different t . 
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Fig. 4. The reflection coefficient for LM 01 -mode at port A 
for different h . 


Fig. 6. Reflection coefficients for ports B, C for different 
septum positions. 


The results of comparison of the scattering parameters in 
the cases of symmetrical (solid line) and asymmetrical (dash 
line) septum position in the interaction region are presented in 
Figs. 5, 6. The step dimensions were chosen as h = 0,5b , 

t-0, 2b h . 

Correctness of the numerical simulation was monitored by 
checking the balance of power in the ports of the T-junction. 
Accuracy of 10‘ 8 was achieved in the considered frequency 
band using 12 modes. It takes less than a second to compute 
100 frequency points using the proposed method. The 
accuracy of the results those obtained by Ansoft HFSS is of 
the order 10“ 2 ...10” 3 and the computation time is of order 
30-50 s. The algorithm allocates memory for about 20 
complex matrices of size 12x12. 

The results obtained by our method are compared with the 
data those obtained by Ansoft HFSS. Agreement within the 
bounds of HFSS method accuracy is achieved for all the 
considered cases. 



4 6 8 

f, GHz 

Fig. 5. Transition coefficients into ports B, C for different 
septum positions 


IV. CONCLUSIONS 

A universal method of rigorous analysis of metallic 
inclusions of arbitrary cross-section located in the interaction 
region of waveguide T-junction is proposed. This model can 
be used to optimize scattering properties of wave guide 
devices using rods or septums placed in waveguide junctions. 
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